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We show that an Anderson Hamiltonian describing a quantum dot connected to multiple leads 
is integrable. A general expression for the non-linear conductance is obtained by combining the 
Bethe ansatz exact solution with Landauer-Biittiker theory. In the Kondo regime, a closed form 
expression is given for the matrix conductance at zero temperature and when all the leads are close 
to the symmetric point. A bias-induced splitting of the Kondo resonance is possible for three or 
more leads. Specifically, for A'^ leads, with each at a different chemical potential, there can he N — 1 
Kondo peaks in the conductance. 
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Introduction. Since the first prediction [l| and real- 
ization of Kondo physics in a quantum dot (QD) 0, 
nonequilibrium effects on the Kondo resonance due to a 
finite bias voltage across the dot have attracted increas- 
ing attention. In the experiments, the zero bias peak of 
the differential conductances has been observed as a sig- 
nature of the Kondo effect on electron transport through 
a QD. In the unitary scattering limit, observations of 
perfect transmission 0, Q provide further evidence for 
the Kondo effect in QDs. The nonequilibrium density of 
states (DOS) of the dot has been predicted ^ to exhibit 
a splitting of the Kondo peak due to a bias voltage ap- 
plied between the source and the drain. This splitting 
has not been observed in transport measurements. To 
observe the splitting of the Kondo resonance by a finite 
voltage bias, an experiment with extra leads (a, has 
been proposed. Very recently, such a splitting was ob- 
served in an experiment 8] where a three-lead setup was 
employed. 

In a conventional bulk Kondo system (e.g., a mag- 
netic impurity in a metal), there is a single chemical po- 
tential and the Kondo resonance in the DOS appears at 
the Fermi energy due to the formation of a singlet be- 
tween the local moments of the impurity and the con- 
duction electrons. If the impurity has available a second 
conduction band to form singlet states, a second Kondo 
resonance in the DOS might be expected to occur at the 
chemical potential of the second conduction band. The 
splitting of the Kondo resonance of a QD by the differ- 
ential chemical potentials of the two leads then seems 
to be reasonable. However, it is not still clear why the 
differential conductance has only a single peak at zero 
bias in experiments with two leads. Thus there arises 
a fundamental question associated with a Kondo reso- 
nance in a system with several chemical potentials that 
can be fabricated in nano-scale electronic devices: why 
the split Kondo peaks have not been seen in two-lead 
systems? To help answer this question we consider a QD 
coupled to multiple leads. The QD is described by an 
Anderson model generalized to a multiple-lead one. It 
will be shown that the multiple-lead Anderson model is 



integrable and exactly solvable by a unitary transforma- 
tion and the Bethe ansatz [l^, [ill ■ By using the 
exact solution, a general expression for the conductance 
of the A^-lead system shows that the Kondo resonance at 
equilibrium is split into N — 1 peaks by increasing the 
difference between the chemical potentials of the differ- 
ent leads. This then clearly shows why only a single peak 
of the conductance occurs in the two-lead system. 

Model. We consider an Anderson model in which N 
leads are coupled to the QD, as in Fig. 1. The leads are 
described under the unfolded formalism with fermions. 
Within this formalism, fermions incident on the dot {x = 
0) from a lead {x < 0) are scattered away from the dot to 
leads (x > 0). In the continuum limit, the multiple-lead 
Anderson model Hamiltonian is given by 

H ^ -i ^ dx cl^^{x)dxCma{x) + Sdd^dq 

N 

+Un^n^ + J2 K^(cL(0)d.+ h.c.), (1) 

rn—l:a 

where rio- — d^d^ is the number of electrons of spin a 
on the dot and U is the onsite Coulomb repulsion. Cmo 
and da are the annihilation operators with spin a for 
electrons in the lead ni and the dot. is a sum over 
the multiple leads [m — 1, - ■ ■ ,N). is the energy level 
on the dot. Here the hopping amplitudes between the 
dot and the lead m, Vm, are allowed to be arbitrary. 
Previously, it has been shown that, for the N — 2 
unitary (Bogoliubov) transformation can be used 
to transform the Hamiltonian to a single-lead Anderson 
Hamiltonian Q|. We now generalize this to the case of 
general N. To do this, one performs a unitary trans- 
formation, c = Ujv c, for the lead electrons, where 
c = (ci,-- - ,CAr) and c = (ci,-- - ,cn)- The compo- 
nents of the N X N matrix Uat are a function of the 
hopping amplitudes, Vm- Uat should satisfy UjyUjv — 

Vm, and (u) 

Z^m Ki[UAr]m,„/ = Vr for m' ^ 1 and for m' 7^ 1, one 
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FIG. 1: A quantum dot (QD) coupled to A''-multiple leads. 
Vm is the tunneling amplitude between the m-th lead and the 
QD. ^Lm is the chemical potential of the m-th lead. The leads 
are presented under the unfolded formalism. 



obtains the one-lead Anderson Hamiltonian and — 1 
free ferniion Hamiltonians. Then a, N x N unitary ma- 
trix for the multiple leads has a form satisfying with 
[UatIi™ = Kn/r and r = Em^^- For > 2, ac- 
tually, there are more freedoms to choose a unitary ma- 
trix. The freedoms give us different matrices for a unitary 
transformation acting only on (c2, • • • ,cn), but leaving 
Ci invariant, which does not affect the physics. 

As a consequence, the unitary transformation satisfy- 
ing such conditions decomposes the multiple-lead Hamil- 
tonian into A^ independent sub-Hamiltonians, Hm, as 



H 



where 



'J — OC 

+Un^ni + Vf {c\^{0)d^ +h.c.) ,(3) 

i^m= / cl^a{x)da,Cmaix) for TOS [2, A^] . (4) 

This is a generalization of the N ~ 2 case treated in 
Ref. 01 . The transformed Hamiltonian can be solved 
exactly because the sub-Hamiltonian, Hi , is the one-lead 
Anderson model that is exactly solvable via the Bethe 
ansatz [H IH El El . 

Integrable excitations and scattering amplitudes. The 
scattering amplitudes of electronic excitations off the QD 
coupled to the A^ leads can be calculated based on the 
exact solution of Hi. In the transformed A^ leads, the 
integrable excitations, {tpm}, will scatter off the dot with 
some pure phase shift with spin cr, 5f (e), where in par- 
ticular S^{e) = for m e [2, A^]. With the unfolded 



formalism, the scattering can be described by the rela- 
tion 



^,n{x > 0) = e'^'^^.mix < 0). 



(5) 



Equation JSJ leads to the scattering amplitudes 5'^„„/ (e) 
of electronic excitations, {ipm}, of energy e between 
leads in the multiple-lead system. Assuming the rela- 
tion ipm = J2mm'i^ N]m,m'4'm' , the Scattering matrix is 
straightforwardly given by 

2 



'S'mm'(e) = ^rnm' + 2irwe* 2 sin -J- , 



(6) 



where Tmrn' — [UPiU~^]mm' and P is a polarization 
matrix: [Pm]rmn = 1 and other entries are zero. For 
m ^ m' , S','^„j/ is a transmission amplitude T^^, from 
to' to TO. For m — m' , corresponds to a reflection 

amplitude from to to to. From Tmrn' = ^m'm, 

^mm'(s) = ^m'm(^) automatically preserved. 

Differential matrix conductance. The current and the 
conductance through the QD can be obtained by the 
Landauer-Biittiker theory 15] for quantum transport 
through nano-devices. To describe scattering away from 
the Fermi energy and calculate the differential conduc- 
tance, we employ an ansatz E3 verified in Refs. [T^ll^. 
The ansatz allows us to use the in-equilibrium scatter- 
ing matrices to calculate the contribution to the current 
of any given excitation. Konik and coworkers discussed 
the details of the implementation of the nonequilibrium 
computation in Ref. T^. With T^^,{e) = T^,,^{e), at 
zero temperature, the current in lead to is given by 



ra' ^m-.a 



de 



(7) 



(2) where jim is the chemical potential at the lead to and 



= 4r^ 



(8) 



To determine Si, we solve Hi via the Bethe ansatz 
for the one-lead Anderson model. The integrability of 
Hi leads to a set of quantization conditions identical to 
that of the one-lead Anderson model. Single particle ex- 
citations with momenta {kj} are identified by an appro- 
priate basis. Scattered particle eigenstates from the dot 
picks up the bare phase S{k) = — 2 tan~^[r/(A: — Sd)]. 
Calculating two particle eigenstates makes it possible to 
get the scattering matrices of excitations. The scattering 
matrices satisfying a Yang-Baxter relationship are iden- 
tical to that of the one-lead Anderson model. Then a 
set of A^e multi-particle eigenstates carrying total spin 
Sz = — M should satisfy the quantization condi- 

tions [la El El as 



rjikj L-\-iS{kj) 



M 

n 



ff(fcj) - Ag -I- i/2 
\ 5(%) - Aa - i/2' 



M 



1 r Aq, — A;3 -I- z 
;!" Aa — A/3 — i 

/3=1 



n 



g(fcj) - Xg- i/2 
\ 9i.kj) - Aq, -I- i/2' 



(9) 
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where g{k) = {k - £d - U/2Y/2UT and M character- 
izes the spin projection of the system with the auxihary 
parameters, {Aq}. For Ed > then, N^. total mo- 

menta fc's form an particle ground state configuration. 
Ne — 2M of iVe momenta fc's is real and 2M is complex 
via M real Aq's. The 2M complex momenta are give n by 
fc± = x(Aa) ± iy{K) with x(A) = U/2 + Sd - VUT[X + 
(A2 + l/4)i/2]i/2 andy(A) = \/Z7r[-A+(A2 + l/4)i/2]i/2. 

According to Andrei's procedure for determining the 
momentum, p, of an added electron in a periodic system 
of size L 19], the quantization condition of the system 
leads to p — 27rn/L. Contributions to the momentum 
come from the bulk of the system and the dot: 

p = 2Tin/L = pbuik + Pd/L. 

The dot contribution scaled by the size of the system 
is identified with the scattering phase of the excitation 
off the dot, which gives the relation between the phase 
and the momentum from the dot as 5i = Pd- In adding 
an electron with spin a to the system, then, the elec- 
tron scattering phase shift has two contributions from 
the charge, p^ ^ and the spin sectors, p^ , as given by 



(10) 



The electronic scattering phase shifts are related to the 
density of states Pd(k) and (Td{\) by the equations: 

Pl{k) = ^(fc)+ /'dA[0i(.9(fc)-A)-27r]arf(A), (11) 



(fc) = ^(fc)+ / dA'[02(.g(fc)-A')-27r]ad(A') 



+ / dfc[0i(A-5(fc))-2^] pd(fc), (12) 



where 5 = 2Re[5(a;(A) -I- iy{\)\- q/B are the Fermi 
surfaces of the seas of fc and A excitations while q is 
related to the band cutoff, D. Here 0i_2 for describ- 
ing the dot momentum should be chosen to ensure that 
Pd^k — > — oo) = pf(A — > oo) = 0. Moreover, the dot 
momenta are simply related to the dot density of states: 

dkp'jik) = 2npd{k), and ^apI (A) = -2nad{\). (13) 

Integrating the density of states gives us the dot mo- 
menta. Consequently, the scattering phase shift is given 
by 

/B 
dkpd{k) + 2TT / dX'adiX'). (14) 
-D Jq 

This phase shift satisfies the Langreth-Friedel sum rule, 
6i = 2TTna, relating the phase shift to the total number 
of electrons Ud in the dot |23| ■ 

To obtain the matrix conductance of the multiple- 
lead system away from the symmetric point (e^ — pm = 



— U/2), we need to do a numerical calculation for the as- 
sociated integral equations. But at the symmetric point 
the scattering phase shift is obtained by using an ex- 
act expression for pd{k < 0) and a direct relation 
between the phase shifts for the electron with spin —a 
and the hole with spin a from a property of electron-hole 
transformation based on the SU(2) spin symmetry The 
phase shift is given by |3 



(5i(e) = -TT— sin ^ 



4Ti,„-7r^(e-Mrn)^ 
4Tl,„ + 7r2(e-^„0' 



C(e), (15) 



where the Kondo temperature for a lead at chemical po- 
tential pm is 



exp 



2rf7 



[iSd-Pni){ed-p,n + U)-r^ 



Here, C(e) does not give any significant phase shift 
when the Kondo energy scale is much smaller than the 
Coulomb interaction U. For \prn — Pm'\ U, we can 
assume all of the leads are at the symmetric point. This 
makes it possible to take into account the essence of 
the physics associated with the splitting of the Kondo 
resonance in a multiple-lead system. Then one can 
obtain a simple expression for the matrix conductance 
{G,nni' ^ -ed^^^Jm) from Eq. 0, (|H1) and The 
matrix conductance in the multiple-lead Kondo-dot sys- 
tem is given by 



(16) 



G =-4Gor: 



2 

mm' 



^2 
1 + ^ 



Pr, 



Pn 



4 VT, 



Jl -1 



,(17) 



where Gq — 2e'^ /h is the quantum of conductance, and 
Fmm' = KnKn'/r- This multiple- lead matrix conduc- 
tance is the generalized expression of the conductance 
for the two-lead Kondo-dot system. It reduces to the 
conductance in the two-lead system Ts^. For a sym- 
metric coupling {Vi — ■ ■ ■ — Vjv) and pi = ■ ■ ■ — pN, 
the matrix conductance is Gmm/Go = 4:{N — 1)/N'^ and 
Gmm'/Go — —{2/NY. The resultant matrix conduc- 
tance agrees with that of a multi-lead quantum point- 
contact for free fermions j2]| . This unitary scattering 
limit shows the Fermi liquid nature of the multiple-lead 
Kondo-dot system. 

Note that the multiple-lead matrix conductance in Eq. 
(|16|l and (|17|l shows clearly that a conductance peak for 
the transmission from m to m! is developed when the two 
chemical potentials are tuned to be equal, pm = Pm'- 
As the chemical potential difference increases, the am- 
plitude of the conductance decreases. In a TV-lead sys- 
tem, if every chemical potential has a different value, the 
conductance Gmm versus pm has a total of the TV — 1 
conductance peaks, one at each of the other chemical po- 
tentials. The amplitude of the conductance Gmm' versus 
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fj,m, has its maximum value for /Zm = fJ-m'- The maxi- 
mum values of Gmm'^s have a one-to-one correspondence 
to the conductance peaks of Gmm- This behavior of the 
conductances implies that electrons from each lead par- 
ticipate in screening the local moment of the dot and 
take part in forming a single Kondo resonance at equi- 
librium. Increasing the difference between the chemical 
potentials, the electrons from each of the N leads have 
their own Kondo resonances with the dot. Each reso- 
nance is characterized by a Kondo temperature, T^.m, 
depending on the value of the chemical potential of the 
lead. Since each lead creates a single lead-dot Kondo 
resonance, the A^-lead system has N lead-dot Kondo res- 
onances. If the chemical potentials of two of the leads 
are adjusted to be equal then the two Kondo resonances 
corresponding to these leads merge together in Gmm'- 
Then this results in only a single transmission peak in 
the conductance Gmm- Therefore, an electron transport 
measurement in the two-lead system is able to capture 
only the single transmission peak even though there are 
two lead-dot Kondo resonances created by the two leads. 
Hence, the two-lead system is not a good probe to observe 
the splitting of the Kondo resonance by finite biases. 

Three-lead and four-lead system. Before proceeding to 
the conclusion, we discuss the conductance for the three 
leads [N = 3) and the four leads {N = 4). The unitary 
transformation for the three-lead system is given by the 
unitary matrix; 



U, = 




Vi V2 V3 

V2 a b 
V3 b c 



(18) 



+ viVr)/j, b = {-V1V2V3 



where a = 

V2V3VT)/-/, and c = i-ViVi ViVr)/-/ with 7 = 
V2 +V^. It can be obtained explicitly under the neces- 
sary condition we discussed above. Similarly, the unitary 
matrix U4 for four leads can be determined. 

We plot the conductance G33 as a function of /13 for 
A'^ = 3 and the conductance Ga as a function of /Z4 for 

= 4 in Fig. 121(a) and (b), respectively. When all the 
leads are at the same chemical potential (A/i — 0), the 
amplitude of the conductance is shown to be reduced as 
the number of leads increases. The maximum amplitudes 
are G33/G0 = 8/9 and G44/G0 = 3/4. As the differ- 
ence between the other chemical potentials, A/z, become 
larger than the Kondo temperature T]^ at equilibrium, 
the single peak at A/z = splits progressively into two 
and three peaks for three and four leads, respectively. 
Figure|21(a) shows that for A/x ~ 2T^, the amplitudes of 
the split peaks reduce to around half the value of that of 
the equilibrium Kondo peak (A/i = 0). The suppression 
of the Kondo resonance is on a voltage scale T^. This 
behavior agree qualitatively, but not quantitatively, with 



the experimental results in Ref. j^]- 

Summary. By using a unitary transformation and the 
Bethe ansatz, the multiple-lead Anderson model is shown 
to be integrable. A general expression for the matrix con- 
ductance from the integrability has been obtained. The 
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FIG. 2: Splitting of the Kondo resonance by multiple leads, 
(a) Conductance G33 as a function of the chemical poten- 
tial /13 for a quantum dot symmetrically coupled (i.e., Vi = 
V2 = V3) to the three leads (A'^ = 3), A/i — /i2 — /ii. (b) 
Conductance G44 as a function of the chemical potential /i4 
for a quantum dot symmetrically coupled to the four leads 
(A'^ = 4), Afi = /i2 — /ii = /ia — /i2. The temperature is zero, 
U = 100. or, and £d = ~3.0r. is the Kondo temperature 
at equilibrium. Here, at equilibrium, all chemical potentials 
are set to be zero. 



conductance for the N-\ead system shows N — 1 split 
Kondo peaks located at — 1 different chemical poten- 
tials. This shows that a Kondo-dot system with multiple 
leads provides a good probe to observe the nonequilib- 
rium effects on the Kondo resonance by a voltage bias in 
transport measurement. 
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